arXiv:l506.00544vl [cond-mat.stat-mech] lJun2015 


Closing the hierarchy for non-Markovian magnetization dynamics 


J. Tranchida a,b , P. Thibaudeau a , S. Nicolis b 
°CEA/DAM/Le Ripault, BP 16, F-37260, Monts, France 

b CNRS-Laboratoire de Mathematiques et Physique Theorique (UMR 7350), Federation de Recherche ’’Denis Poisson” 
(FR2964), Departement de Physique, Universite de Tours, Parc de Grandmont, F-37200, Tours, France 


Abstract 

We propose a stochastic approach for the description of the time evolution of the magnetization of nano¬ 
magnets, that interpolates between the Landau Lifshitz-Gilbert and the Landau-Lifshitz-Bloch approxi¬ 
mations, by varying the strength of the noise. In addition, we take into account the autocorrelation time of 
the noise and explore the consequences, when it is finite, on the scale of the response of the magnetization, 
i.e. when it may be described as colored, rather than white, noise and non-Markovian features become 
relevant. We close the hierarchy for the moments of the magnetization, by introducing a suitable truncation 
scheme, whose validity is tested by direct numerical solution of the moment equations and compared to 
the average deduced from a numerical solution of the corresponding stochastic Langevin equation. In this 
way we establish a general framework, that allows both coarse-graining simulations and faster calculations 
beyond the truncation approximation used here. 


1. Introduction 

At atomistic length-scales, relaxation processes 
towards equilibrium of the magnetization in mag¬ 
netic systems rely on complex interactions between 
spin, electron and lattice subsystems [El, 9- The 
time-scales for these processes are extremely short 
(t « 100 ± 80ps for a ferromagnetic gadolinium 
system 0 ) , so a first approach for their description 
was to perform stochastic simulations of interacting 
spins while neglecting the memory effects for the 
noise (Markov’s hypothesis) [b, fj. However, recent 
experimental breakthroughs are pushing the time- 
scales that can be probed for the magnetization dy¬ 
namics to pico- and even femtosecond resolution[7|. 
Furthermore, a self-consistent approach for the 
short-range exchange interaction, that is responsi¬ 
ble for local spin alignment, indicates that it would 
be equivalent to a very high value of the local mag¬ 
netic field (more than 80 MAm -1 ). To take into ac¬ 
count this contribution in simulations, it seems to 
impose time-steps at least one order of magnitude 
shorter than all the physical time scales (dynamics, 
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relaxation, etc.) of the system. Under these condi¬ 
tions, numerical simulations of the stochastic effects 
in such magnetic systems are raising the question 
of the validity of the Markovian assumption that is 
commonly used, that is to have short memory in 
the sense that the correlation time is very short [g]. 
In order to test this assumption, by taking mem¬ 
ory effects explicitly into account, a colored-form 
for the noise has to be considered. In spin systems, 
the noise is represented by uj as a random vector 
whose the components, w, are Gaussian random 
variables with zero mean and a the finite correla¬ 
tion time, which is encoded in its two-point function 
as follows: 

= fd, exp (22——^ (1) 

where D is the amplitude of the noise, and 
r is the correlation time. We remark that 
lim y exp = 5{t — t') which allows the 

white-noise expression for the correlator to be re¬ 
covered in this limit. 

To compute the moments of the magnetization, 
subject to such noise, we adapt the formalism de¬ 
veloped by Shapiro and Loginov 9 to the magneti¬ 
zation dynamics. A new hierarchy of equations for 
the moments is then deduced and appropriate clo- 
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sure relations allow us to solve them directly. The 
solution is compared to numerical simulations of 
the corresponding stochastic Langevin equations, 
that are performed in the vicinity of the white noise 
limit. 


2. Moment equations and closure of the hi¬ 
erarchy 

The Einstein summation convention is adopted, 
with Latin indices standing for vector components, 
eijk is the anti-symmetric Levi-Civita pseudo¬ 
tensor and Si are the components of the normal¬ 
ized magnetization vector (|s| = 1). The stochas¬ 
tic Landau-Lifshitz-Gilbert (sLLG) equation can be 
written as follows: 


that have to be derived. By applying the Shapiro- 
Loginov method [l2j to the nine components of 
the ( uiiSj ) matrix at the same time, we obtain the 
equation: 


dt 



1 

T 


(tiiSj) 


(5) 


Injecting the right-hand side of Eq.© in Eq.©, 
equations for the second-order moments (tiiSj) are 
obtained. 

The equations for the moments (si.Sj) are now 
required. Assuming that the following identity 


djsiSj) 

dt 



( 6 ) 


dSi 1 r i - 

~Q^ - l + \2 e ijkSk [Uj + u>j 

^ Td^)s m )] (2) 

Here w is the effective precession frequency of the 
magnetization, that is assumed to be constant, A 
is the damping coefficient, and w represents the 
stochastic noise contribution, whose components 
are drawn from a Gaussian distribution with zero 
mean and a colored form for the second order corre¬ 
lator eq. ©. This is the fundamental equation for 
the dynamics of a small magnet in contact with a 
stochastic bath, which we shall take to be thermal 
henceforth. It has the formof a Langevin equation, 
with multiplicative noise sa¬ 
lt can be simplified by keeping only one random 
torque on the right-hand side: 


dsj 

dt 


y ^2 \pJj + Cjj Aej/ m CU/S m ] (3) 


Indeed, up to a renormalization of the noise, it can 
be shown that this equation generates a stochastic 
dynamics that is equivalent to that of eq. ®.H]- 
These microscopic degrees of freedom give rise 
to average quantities, that probe statistical ensem¬ 
ble behavior. These can be operationally defined 
by taking statistical averages over the noise 0 of 
eq. ©• Denoting by (.) this average, we have: 


d(si) 

dt 


1 ^2 [ e ijkWj(sk) + CijkfajSk) 

A eijkCjlm UJl(s k S m >] (4) 


This average magnetization dynamics relies on 
higher order correlation functions of noise and spin 


holds, the time derivatives of dsi/dt are replaced 
by the right-hand side of Eq. © and so forth. 

This leads to an additional set of nine equations 
for the second-order moments of the microscopic 
spin degrees of freedom, that are closely related to 
those written by Garanin et al. (see Eq.(6) in ref. 
jl3l p and Gracia-Palacios et al. (see Eq.(2.10) in 
ref. 0). 

The complete system of equations, therefore, 
takes the following form: 


d{si) 

dt 

d{tiiSj) 

dt 
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dt 
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Imn^m (SiSiS n }) + (i j)( 9 ) 


Due to the non-linearity of the sLLG equation, an 
infinite hierarchy arises 15. these equations for 
the second-order moments depend on third-order 
moments, and so on. In order to solve the system, 
closure relations need to be found, that express the 
third-order moments in terms of the second order 
and first order moments. These relations can be 
deduced from the definition of the corresponding 
cumulant moments. 

The double bracket notation ((.)) stands for the 
cumulant of the stochastic variables H3 , and for 
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any stochastic vector x, one has: 

(( XiXjXk)) = ( XiXjXk) - (Xi)(XjX k } 

-( Xj)(XiX k } - ( X k )(XiXj) 

+2 (xi)(xj)(x k ) (10) 

Assuming that the third-order cumulants vanish for 
each stochastic variable (meaning that the corre¬ 
sponding distribution is assumed Gaussian), and 
applying the relation m, the following closure re¬ 
lations are deduced: 


(' UiSjSk) = (Sj)(u>iSk) + (Sk)(u>iSj) (11) 
(tiiCjkSi) = —Sik(si) (12) 

T 

(SiSjSk ) = { s i){ s j s k) + ( Sj)(SiSk) 

+ (Sk)(SiSj) - 2(Si)(Sj)(s k ) (13) 


The second closure relation arises from the prop¬ 
erties of the noise for the same time: indeed, one 
has (cDj) = 0 and {oJiOJj) = 5y when (Dj and u>j 
are considered at the same time. Introducing those 
relations in the system ©, we finally have: 


d{si) 
d t 


d{u>iSj) 

dt 


9(siSj) 

dt 
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^ [ejkiUk(&iSi) 

D / \ 
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+ < s n )(uiSi ))] (15) 

1 ^2 i e jklUk(SiSl) 

+Cjki ((si)(uikSi) + (si)(QkSi)) 

^tjkltkmnUm {slS n ) 

+ (sz)(SiS ra ) + ( S n )(SiSl ) 

- 2 (si)(si){s n ))} + (it) j) (16) 


This set of twenty-one differential equations is now 
closed and can be integrated simultaneously to de¬ 
scribe the averaged dynamics of the effective mag¬ 
netic system connected to a bath. It should be 
noted that the last two equations involve both, lon¬ 
gitudinal and transverse, forms of damping. In con¬ 
sequence thereof, the first equation, that describes 
the evolution of the averaged magnetization, has a 
Landau-Lifshitz-Bloch form, with both, transverse 
and longitudinal, damping terms. 

In the following sections, this model will 
be called the dynamical-Landau-Lifshitz-Bloch 
(dLLB) model. 


3. Numerical experiments 

Numerical experiments are performed in order to 
test the consistency of our closure hypothesis. 

We first focus on the limiting case of small cor¬ 
relation time, i.e. the vicinity of the white-noise 
limit. Indeed, if r is very short (the same order 
of magnitude as our integration time), the associ¬ 
ated relaxation can no longer be resolved, and loses 
its physical meaning. Under those conditions, the 
dynamics is Markovian, and our dLLB model is ex¬ 
pected to match the results obtained from an aver¬ 
aged Markovian stochastic magnetization dynam¬ 
ics. 

In order to realize these comparisons, a reference 
model is built from the sLLG (Eq. © ) with a white- 
noise random vector. One thousand realizations of 
this stochastic equation with the same initial con¬ 
ditions are generated, and a statistical average is 
taken over them. The dLLG system is solved by us¬ 
ing an eight-order Runge-Kutta integration scheme 
with adaptive steps constructed by the Jacobian of 
the whole system. 

Fig. (© shows the average, obtained from the 
sLLG equation, and the dLLB model for the same 
value of the noise amplitude D. It appears that the 



Figure 1: Relaxation of the averaged magnetization dy¬ 
namics about an external field in the ^-direction and in 
the case of small correlation time for the noise. The solid 
lines plot the stochastic case (1000 averaged repetitions), 
dots are for the dLLB model. Simulations parameters: 
{D = 1.6.10 2 ; A = 0.1; r = 10 3 ; 63 = (0, 0,7r/5) }. Initial 
conditions :s(0) = (1,0,0). 

magnetization dynamics obtained with our dLLB 
model fits very well the averaged sLLG results in 
the case of small correlation time for the noise. 
Moreover, an essential feature of the Fig. © is 
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that our model does behave like an LLB equation 
and can describe the variation of the magnetiza¬ 
tion norm, that does not remain constant. This is 
due to the longitudinal damping component that is 
generated by the average effects in presence of the 
noise. 

The second order moments, also, play a major 
role in the evolution of the averaged magnetization 
dynamic equations, and the diagonal components 
of the tensor (sjSj) are plotted in Fig. ©• 

Good agreement is also found between the two 



Figure 2: Relaxation of the diagonal components of the cor¬ 
relation tensor ( S{Sj ). Again,the solid lines plot the stochas¬ 
tic case (1000 averaged repetitions) and the dots are for 
the dLLB model. Initial conditions: (si(0) 2 ) = 1, and 
(*(0)^(0)) = 0V(i,j)^ (1.1) . 

models for the diagonal components of the correla¬ 
tion tensor ( SiSj ). Those results allows us to test 
the validity of the equality Eq. ©, and the closure 
relations Eqs. (111112113b . 

4. Conclusion 

Statistical average and the application of the 
Shapiro-Loginov method allow us to derive a new 
formalism for the magnetization dynamics when 
a colored form for the noise is appropriate. Our 
model has proven to be effective for the description 
of an average magnetization dynamics when small 
values of the time-correlation r are considered. 
In this particular case, the results of Markovian 
dynamics are recovered and very good agreement 
was found between our dLLB set of equation and 
an averaged set of LLG equation with a white-noise 
process. Moreover, this work demonstrates that 
the application of statistical average to a magnetic 


system in contact with a thermal bath leads to 
the direct appearance of a longitudinal form of 
damping, and to equations consistent with the 
Landau-Lifshitz-Gilbert model. In future work, 
the effect of longer correlation times for the noise 
will be studied in detail, and comparison between 
our dLLB models and atomistic magnetization 
dynamics in more complex situations (exchange 
interaction, anisotropy, superparamagnetism, etc.) 
will be presented. 
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